Confinement and lattice QED electric flux-tubes simulated with ultracold atoms 
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We propose a method for simulating 2+1-d compact lattice quantum-electrodynamics (QED), 
using ultracold atoms in optical lattices. In our model local Bose-Einstein condensates' phases 
correspond to the electromagnetic vector-potential, and the local number operators represent the 
conjugate electric field. The well-known gauge invariant Kogut-Susskind Hamiltonian is obtained 
as an effective low energy theory. The field is then coupled to external static charges. We show that 
in the strong coupling limit this gives rise to 'electric flux-tubes' and to confinement. This can be 
observed by measuring the local density deviations of the BECs, and is expected to hold even, to 
some extent, outside the perturbative calculable regime. 



Free quarks are not found in nature; This is due to the 
mechanism of confinement A lot of theoretical progress 
in this area has been achieved - either in the lattice Eu- 
clidean approach by Wilson [1 j, in nonperturbative meth- 
ods by Polyakov [2] or using the lattice Hamiltonian for- 
malism, by Kogut and Susskind [3l 0] . 

Although gauge theories can be latticized either in a 
compact (nonlinear) or noncompact (linear) manner, the 
compactness is essential to the confinement mechanism 
[5]. It has been shown that in an abelian 3+1 compact 
lattice gauge theory, a phase transition is supposed to 
take place between two phases - the Coulomb phase for 
small couplings, which exhibits the 'regular' V (R) oc 1/R 
static potential between two ^-separated static charges, 
and the confining phase, for which the static potential is 
linear in the distance between the charges - V (R) oc i?, 
for large values of the coupling constant [4 j. (Nonabelian 
theories, on the other hand, confine for all values of the 
coupling constant). However, for an abelian 2+1 com- 
pact lattice gauge theory, confinement was shown to take 
place for all the values of the coupling constant, due to 
nonperturbative effects of instantons [2j [SHT] . Consid- 
ering thermal effects as well, even in 2+1 dimensions a 
phase transition to a Coulomb phase exists for T > 

Hi. 

The mechanism responsible for confinement is believed 
to produce an 'electric flux tube', connecting two static 
charges in the confining phase, which is hard to measure 
directly. It requires measuring the force and/or potential 
between two static charges. If one wishes to observe the 
phase transition, the coupling constant has to be varied, 
which poses another difficulty. A quantum simulation of 
such a model could allow a direct test of the confinement 
mechanism and the phase transitions. 

Quantum gases of ultracold atoms, implemented in op- 
tical lattices [10] , provide models with highly controllable 
parameters and offer a natural playground for the simu- 
lations of such models. Quantum simulation approaches 
of various kinds and aspects of compact U(l) pure gauge 
theory, in cold gases and other systems, have been pro- 
posed by several authors: In [IT] , an effective theory of 



U(l) spin liquid in pyrochlore was discussed; in [12], using 
a molecular state in optical lattices, an effective theory 
of ring exchange was derived, and it is, in the limit of no 
hopping, a U(l) lattice gauge theory, with a Coulomb 
phase; in [13], emergence of 'artificial photons' and a 
Coulomb phase in an effective theory based on dipolar 
bosons in an optical lattice were shown; and in [14], a 
possibility to simulate a spin U(l) pure gauge theory as 
a low energy theory with a system of Rydberg atoms was 
presented. 

In this letter, we suggest a method for simulating com- 
pact QED with cold atoms in optical lattices, which 
should enable a direct observation of 'electric flux tubes' 
that emerge in the mechanism of confinement. In our 
model, the vector potential and its conjugate electric field 
are represented by the local condensate phase operators 
and their conjugate number operators. These observables 
'live' on the links of a two- or three-dimensional optical 
lattice, and hence each link of the lattice is here repre- 
sented by a separate Bose-Einstein condensate. In order 
to obtain the QED Hamiltonian, one has to generate cer- 
tain two- and four-body interactions between the conden- 
sates, that manifest local gauge invariance. In order to 
avoid the hopping processes of an ordinary Bose-Hubbard 
model, we introduce a four-species two dimensional set- 
up (fig. 1). The condensates are located on the links 
of a lattice - each species on a different link - and over- 
lap at the lattice's vertices. Hence, condensates of the 
same type are spatially separated, as depicted in figure 
1, causing the attenuation of hopping processes. Next we 
use Raman transitions and two atom scattering processes 
in order to create special 'diagonal' hopping and nonlin- 
ear interactions. As we show, in this new setup a certain 
choice of parameters gives rise to gauge invariance in the 
low energy sector, hence compact QED emerges as an 
effective theory. 

To study the effect of confinement within this setup, we 
can introduce two spatially separated effective 'charges' 
by creating local deformations of the trapping potential 
at the position of the 'charges' at the relevant vertices. 
We then expect that the local atomic densities, within 



2 



the QED parameter regime should manifest the effect 
of confinement by the appearance of a flux-like tube of 
alternating atomic density deviations along the line con- 
necting the 'charges' (fig. 2), while such a 'flux tube' will 
not appear outside the QED parameters regime. Other 



Sij is Kronecker's Delta, gij are the s-wave scattering co- 
efficients and Qij are Rabi frequencies. It contains the 
following parts: The 'free' Hamiltonian of each species: 
Wq (x) = — -f y i ( x ), where Vi (x) is the optical lattice 
trapping potential of the species z; the scattering terms, 
set by the coupling constants g^ (neglecting the three- 
and four-body interactions): (i) self-scattering terms - 
gu = gi, (ii) two-species scattering terms: g±2 = gi\ = 
934 — 943 = 92 along straight lines, and along the diago- 
nals, 6/13 = 931 = gi4 = 941 = 923 = 932 = 924 = 942 93 

(all the other g^-'s are zero); an 'external charges' sim- 
ulating potential, which deforms the lattice potential at 
the vertices and is approximated by a very localized po- 
tential: Vm(x) ]Ta m)n £(x-x m)n ), where a m , n are 

constants (whose value and sign are related to the 'ex- 
ternal charges' and will be determined in the sequel) and 
x m?n is the position of the (m, n) vertex. The laser gen- 
erated 'Rabi terms' ^13 = ^31 = ^14 = ^41 = ^23 = 
^32 = ^24 = ^42 = ^0 couple the condensates to each 
other in a special, 'diagonal' manner, as depicted in fig- 
ure [l] (All the other Q^-'s are zero). Since the minima 
of the same species are far enough apart, the hopping ef- 
fects are solely controlled by the latter 'Rabi terms'. Ex- 
perimentally, our scheme can be implemented by using 
holographic masks techniques [15] in order to generate 
the required optical lattice and using optical Feshbach 
resonances in order to control the coupling strengths g^ 
[T6HT8] ). Raman transitions can be used to control the 
coefficients Qij of the 'Rabi terms'. 

The second quantization wavefunctions of the con- 
densates (taking into account only the lowest band 
excitations) are ^ lj2 = X) a m,n</>m,n (x) , ^> 3 ,4 = 

m,n 

^2 b m,nXm,n (x) where a m , n ,6 m5n are single-mode anni- 

m,n 

hilation operators, annihilating one particle in the ground 
state of the corresponding link (minimum) . Note that be- 
cause of the lattice's structure, not all the values (m, n) 
are included in the wavefunction of each species. We 
assume that the local Wannier functions [10] respect 
the symmetries m>n (x) = (x - x^ 2 n ) , x m , n (x) 
X (x — x^ 4 n ) = (j) (i?x), where R is the appropriate rota- 
tion operator, and that they are chosen to be real [19]. 

Plugging the wavefunctions into the Hamil- 



possible implications of our model will be shortly dis- 
cussed in the summary. 

We bgein with a system of condensates described by 
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the Hamiltonian H = J d 3 x ^ Hij (x), where 

(1) 



r 




FIG. 1. On the left - structure of the lattice. The differ- 
ent condensate species are colored in four colors; the colored 
boxes represent the links (condensates), and there the local- 
ized wavefunctions are concentrated. At the vertex (symbol- 
ized by a cube) the wavefunctions of the neighboring links 
overlap and these are the only overlap integrals which are not 
negligible. On the right - a close-up picture of a single ver- 
tex, showing the various interaction parts of the Hamiltonian 
- scattering and hopping. 



tonian ([I]), one gets, using the above assump- 
tions, that the only non- negligible contribu- 
tions are: 2A + /i = ^ J d 3 x \<p (x — x )| 4 , 

V = J d 3 xr (x - x ) (-^ + Vi (x)) (x - x ), 

V 2 EE Jd 3 x|(/>(X-X )| 2 |(/>(X-X 1 )| 2 , 

V 3 EE /d 3 x|0(x-X O )| 2 |x(x-X 2 )| 2 , 

A m , n ee -^a m ,n |0(x m , n )| and Q 
Q f Q f d 3 xcj)* (x — xo) X (x — X2) (Here the reality of 
the Wannier functions is employed), xq is the position 
of an arbitrary minimum of the potential (due to the 
symmetries), xi is an adjacent minimum in the same 
direction (separated by a single lattice spacing) and 
X2 is an adjacent minimum in the orthogonal direction 
(rotated). In the following we assume that g 2 ,#3 satisfy 
the relation g 2 V 2 = 93^3 — 2 A. 

Let Nj^ n be the local number operators, emanating 
from the vertex (m, n): for horizontal (x) links k = 1 and 
for vertical (y) links k = 2. N T = ^m,n> the total 

m,n,/c 

number of particles, is a constant of motion. We choose a 
subspace by fixing Nt — A/lAT , where Ml is the number 
of links and N > 1). Defining M m>n = 47V + A m>n , 
G m ,n = N^ n + N^ n + Nl_^ n + - M m , n , after 

some algebra, one obtains the Hamiltonian Hq ee H — 



I 

Hij (x) = *J (x) (d l3 (Hi (x) + VW(x)) + (x) + 9 -f*\ (x) *t (x) ^ (x) tf, (x) 
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H R = A E + /" E ( N m,n) 2 - The nearest-neighbor 

m,n m,n,k 

hopping part, which results here from the 'Rabi terms', 
can be written as H R = () ^ (a^n&t^ + am,n&m + i >n + 

ra,n 

Gauss ; 5 Lau>. We wish to obtain a gauge-invariant the- 
ory, and hence we would like to constrain Gauss's law 
on the system. This is satisfied in the 'QED regime': 
A ^> /i and A ^> Qq, in which Hr can be treated 
as a small perturbation. Let us first find the ground 
state of Hq. After expanding the number operator on 
each link around Nq : n = Nq + 5^ n , one ob- 
tains at each vertex: G m , n = S^ n + 5^ n + ^_i jn + 
n-i — ^m,n- Within the subspace of a constant, con- 
served Nt, ^2 5^ n — 0. Neglecting constants of mo- 

m,n,/c 

tion, one can rewrite the Hamiltonian in terms of n : 
H = XE + n E (SLn) 2 = H G + H E . A » M , 



mm,k 



and hence one would like to minimize Hq first. Thus we 
get that in the ground state, the sum of 5^ n 's around 
each vertex equals the A m ^ n of the vertex: This imposes 
a 'modified Gauss's Law' (sum instead of discrete diver- 
gence), and hence the A m)U 's must be integers (positive, 
zero or negative) - this can be set by adjusting the values 
of the a m)n 's in Vm(x). Next, to minimize the entire Hq 
(including He) we would like to choose the lowest 5^ n 's 
which satisfy this constraint. 

Define the sublattices A — {(m,n) : m + n = even}, 
B = {(m,n) : m + n = odd}. Note that for states that 
respect Gauss's Law (for which G min = 0), which 
will later be the physically interesting states, the sum of 
A m , n 's of each sublattice must be zero. This follows from 
adding the G m)n 's of each sublattice, taking into account 
that the total particle number deviation is zero. 

Quantum Rotor Approximation. If we set that at each 
vertex |A m5n | <C Nq, we get that on each link, in the 
ground state of Hq, within our subspace, |(5^ n )| <C 
Nq, and thus, after taking into account the perturba- 
tive corrections, one obtains that on each link a% 

1/2 
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< N . Note that a: 
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mm 



XNo(Nq-I) 
Mo 



(because ^ <C Nq), and hence the two conditions of 
[20] are fulfilled and the Hamiltonian can be approxi- 
mated as a quantum rotor Hamiltonian. Hq,He re- 
main the same, because they are already written in 
the number deviations' notation. Because of the phase- 
number relation of the condensates, 
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n 1 e 3 , , 

±y mmi u m' m' > 

> n and therefore 



we can define phase-only lowering and raising operators, 



1 , a, 



1 



e ^m,n ? which operate on the local number devi- 
ations: a myn \S^ n ) = \& m ^ n - l) etc., and since 
Nq > 1, ^Nq (Nq + 1) ~ _Nq. Thus one gets 

Hr = ( a mm^ln jn + a m,n^ m +i )n + a mm+l^ m m 

mm 

2 m , n +i^|n+i,n + h ' G ')i where ft = ft N . 

Effective Hamiltonian. Let us look again at the eigen- 
states and eigenvalues of Hq- Since [Hq,He] = 0, the 
two Hamiltonians can be mutually diagonlized. The 
eigenstates of He are number states, and we shall use 
this basis to diagonalize Hq as well. Since A ^> \i, Q, 
the Gauss Hamiltonian Hq is much stronger than the 
other two, and therefore one can obtain an effective low- 
energy theory perturbatively [21 j. It is physically rea- 
sonable to derive an effective Hamiltonian by project- 
ing to the ground state manifold of Hq- Let us de- 
note this manifold by M: M = {\M a } : Hq \M a ) 0}. 
One can see that it is the physical subspace of states 
which respect Gauss's law. The perturbative expan- 
sion to second order leads to H e ff = He + Hb, where 



2Q- 
A 



E (ah, 



&m+i,n + h.c. ) is the 



desired gauge invariant four-body plaquette interaction. 

Compact QED analogy. We next relate this model 
to compact QED and discuss the implications. First, 
let us switch to QED-like variables. In order to do so, 
we perform the transformation: E 1 ^ nn = (— l) m+n ^ n , 
Qm,n = (-l) m+n A m , n and 6^ n V(-i) m+n #m,n- Be- 
cause of the transformation of the phases of links ema- 
nating from sublattice B vertices, these links' raising and 
lowering operators have to be swapped. 

This transforms He, which can be identified as 
the 'electric Hamiltonian', to He = /i ^ (E t 

m,n,k 

and H B = ~^E (aL.nSm.n+l&L+l.n^™," + h.C 



^m,n ) 
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Ecos (01+0; 
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m+l,n 



m,n+l 



IS 



the 



magnetic part of the compact QED Hamiltonian (the 
cosine's argument is the discrete curl of m?n , which 
is the magnetic field). Thus we obtained an effective 
low-energy theory whose Hamiltonian is the compact 
QED Hamiltonian, constrained with Gauss's law (which 
is the low-energy constraint): 



= (-l) m+n (divE 



QmmM =0 (2) 



Confinement of external static charges. Define a new 



finite energy scale, L^o 
Hamiltonian to 



4^V 



and rescale the 



H = H e ff/Uo = y ^ { E m,n) C ° S i S m,n + °m+l,n ^m,n+l ^m,n) ( 3 ) 

mm,k mm 
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FIG. 2. An example of the charge and flux configurations, for 
R=2 . The different colors represent the condensate species. 
The upper couple of charges are with QED quantum num- 
bers, and the lower couple with BEC local number deviations 
quantum numbers. Such a 'flux tube' can be embedded, in 
absence of other charges, in a lattice whose other links carry 
E = 0. 



which is the well-known Kogut-Susskind Hamiltonian for 
an abelian lattice gauge theory [2 [4] . From the definition 
of U one gets g 4 = This Hamiltonian has two limits: 

(i) The strong coupling limit: g > 1, or ^ > j. In this 
limit, we can treat Hb as a perturbation to He, and 

(ii) The weak coupling limit: g <C 1, or ^ <C In this 
limit, we can treat He as a perturbation to Hb- In a 3+1 
theory, the strong coupling limit is within the confining 
phase and the weak coupling limit is within the Coulomb 
phase, and a phase transition is expected in between p]- 
[5]. In a 2+1 theory, there is no phase transition and 
confinement is expected to occur for all g > [5j [7] . 

The 'external charges' are limited by the restric- 
tions imposed by the constraint Y A m?n = 

(m,n)eA 

0, ^m,n = 0- subtracting the second constraint 

(m,n)eB 

from the first, one gets ^ Q m ,n = Qtotai = 0. Thus, the 

m,n 

total charge has to be zero. If we add these constraints, 
we get another constraint , Y ^m,n = 0. This constraint 

does not seem to have a QED analogy, but it has to be 
satisfied in our model. 

Consider the case of a system with two unit 'external 
charges', in the strong coupling limit. Thus we seek for 
the ground state configurations of He and treat Hb as a 
perturbation. These 'charges' must be of opposite signs, 
in order to satisfy the charge restrictions. For simplicity, 
we assume that the 'charges' are fixed at the vertices 
(m,n) and (m + R, n). If R = 1, the 'charges' are fixed 
at two vertices of different sublattices, and hence have 
to have the same sign in terms of A min . This, however, 
does not satisfy the charge restrictions, and one has to 
add more 'charges' to the system. This is true, in fact, 
for any odd R. Therefore we shall consider only the case 
of an even i?, for two 'charges' in the system. 

Denote \R) as the state of two such 'external charges', 



zeroth order term is \Rq) = ] [ aj^ n |{0}). This 

mo<m< rn q-\-R 

corresponds to a 'flux tube' from a positive charge to a 
negative one. Thus, in the strong limit, we get, indeed, 
the expected strong coupling linear behavior of the en- 
ergy, 



E(R) 



E{R) = 9 -R + 0(g- 



U 



(4) 



The effect can then be observed by measuring the local 
density deviations 5^ n , which are expected, in the lead- 
ing order, to have a magnitude 1 and alternating signs 
between the two 'charges'. An example for R = 2 can 
be seen in figure [2j When R is too large, the energy of 
Hq is smaller than the energy of such a 'flux tube', and 
then the low energy theory breaks, and a 'flux tube' is 
no longer the state of minimal energy. The low energy 
picture holds as long as the 'flux tube' length satisfies 
R < A//i. 

Outside the strong coupling regime, such perturbative 
calculations are no longer valid. However, in 2+1 dimen- 
sions the confinement holds for all values of g. Hence, the 
effect should be seen experimentally even slightly outside 
the strong coupling limit, although not in the weak limit 
(in order to fit with the Quantum Rotor approximation) . 

It may also be possible to experimentally observe the 
effect of a finite temperature, T > 0, on the model, in- 
cluding a phase transition. 

Extensions of the model In this paper, we have shown 
a method to simulate compact QED using BECs in op- 
tical lattices, as a way to observe charge confinement. 
The suggested model can be extended in several ways: 
More realistic, nonperiodic boundary conditions can be 
imposed (e.g., no 'charges' on the boundary); Using ad- 
ditional condensates (of new species), a 3+1 simulations 
could be achieved. Interestingly, a dynamical charge 
which is minimally coupled to the field can be imple- 
mented using another condensate species. This is equiva- 
lent to a special case of the model of Fradkin and Shenker 
[22] . in which a Higgs field with a 'frozen' radius is cou- 
pled to a U(l) gauge field. [23] . 

To Conclude, we hope that this model could serve as 
one of the building blocks of the bridge into the world of 
dynamic quantum gauge field theories simulations. 
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